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; Abstract 

The possible range of n — n' mixing angle is determined from the transition form factors F m (Q 2 ) 
and F r] r 1 (Q 2 ) with the help of the present experimental data. For such purpose, the quark- flavor 

o 

■ light-cone pQCD framework, where the transverse momentum corrections and the contributions 



9" 
6. 



mixing scheme is adopted and the pseudoscalar transition form factors are calculated under the 



\Q . beyond the leading Fock state have been carefully taken into consideration. We construct a 

o 

phenomenological expression to estimate the contributions to the form factors beyond the leading 



Fock state based on their asymptotic behavior at Q 2 — > and Q 2 — > 00. By taking the quark-flavor 
mixing scheme, our results lead to <p = 38.0° ± 1.0° ± 2.0°, where the first error coming from 



experimental uncertainty and the second error coming from the uncertainties of the wavefunction 
parameters. The possible intrinsic charm component in n and n' is discussed and our present 
analysis also disfavors a large portion of intrinsic charm component in r\ and 77', e.g. \f^, \ < 50 MeV. 
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I. INTRODUCTION 



The light-cone (LC) formalism [jll^S] provides a convenient framework for the relativistic 
description of hadrons in terms of quark and gluon degrees of freedom and for the application 
of pQCD to exclusive processes. Among them, the 77 — 7/ mixing is a subject of considerable 



interest, which has been examined in many investigations, see e.g. 

y 



llltClCOL, W111UL1 lid , 

y, y, 0, y, y, q, 



11 



12 



13] and theoretical ones 
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13, 



experimental ones 
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Some experiments have been done recently, e.g. the new KLOE value of — r ^~* v ^ 



(4.9 ± 0.1 stat ± 0.2 syst ) ■ 10" 3 |s| leads to <p = 41.2° ± 1.2° [l2j; the BES collaboration has 



r[0-7?7] 



announced a new measured value for Rj/q, i.e. Rj/y = = 

4.94 ± 0.40 [13j that 

leads to = 38.8° ± 1.2°. Furthermore, the newly measurements of the form factors F m (Q 2 ) 
and F^ 7 (<5 2 ) at asymptotic region by BaBar collaboration [lOj, Q 2 F v ^(Q 2 )\q2 =U2 Gev 2 = 
0.229 ± 0.030 ± 0.008 GeV and Q 2 F v ^(Q 2 )\ Q 2 =ll2 GeV 2 = 0.251 ± 0.019 ± 0.008 GeV, will 
provide further constraints to the theoretical predictions. 

The pseudoscalar transition form factors F m (Q 2 ) and F r) i 1 {Q 2 ) provide a good platform 
to study the 77 and 7/ mixing effects, which have already been studied in literature by several 



groups 
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21 
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251 ] . However in these calculations, either only the leading Fock- 



state (e.g. Ref.[24j) or only the asymptotic behavior of the form factors (e.g. Ref. [25]) have 
been taken into consideration to determine the mixing angle. As has been pointed out in 
Refs. 



21 



23] , the mixing angle can not be reliably determined without a proper considering 



of the contributions from the non-valence quark states due to the fact that even though 
the higher Fock states' contributions are power suppressed in large Q 2 region, they will 
give sizable contributions to small and intermediate regions. In fact, it has been pointed 



out in Ref. 



26] that the leading Fock state contributes to F p ^(Q 2 )\q2 =0 (P stands for the 



pseudoscalar mesons) only half and the remaining half should be come from the higher Fock 
states. The higher Fock states' contributions in small Q 2 region can not be calculated by 
the perturbative QCD approach due to its nonperturbative feature. Recently, Ref. {2^] had 
constructed a phenomenological expression for the pion-photon transition form factor to 
estimate the contributions beyond the leading Fock state based on its asymptotic behavior 
at Q 2 — > and Q 2 — > 00 (Q 2 stands for the momentum transfer in the process). The 
predicted results for F 7T ^ / (Q 2 ) there agree well with the experimental data in the whole Q 2 
region. In the present paper, we will adopt this newly developed method to estimate the 
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higher Fock states' contributions of the form factors F V1 (Q 2 ) and F r? / 7 (Q 2 ), and then to 
derive a possible range for the mixing angle by comparing the predicted results with the 
experimental data. 

As for the 77 — 77' mixing, two mixing scheme are adopted in the literature, i.e. the 
octet-singlet mixing scheme and the quark-flavor mixing scheme. These two schemes can 
be related with a proper rotation of an ideal mixing angle (9 id = — arctany^ — —54.7°) 
rsl . H ]. A dramatic simplification can be achieved by adopting the quark-flavor mixing 
scheme, especially, the decay constants in the quark-flavor basis simply follow the pattern 
of state mixing due to the OZI-rule 181 ] . Furthermore, by adopting the quark-flavor mixing 
scheme and also by carefully dealing with the higher Fock states' contributions, a naive 
discussion (at the end of Sec. II. C) shows that the value of Q 2 F V1 (Q 2 ) decreases, while the 
value of Q 2 F T] i 1 (Q 2 ) increases, with the increment of the mixing angle 0, so a possible range 
for can be derived by comparing with the experimental data on the form factors F rri {Q 2 ) 
and Frfy(Q 2 Y And then by adopting the relation between the two schemes as shown in 



Refs. 



18 



281 ]. all the three mixing angles 9p, 9\ and #8 involved in the octet-singlet scheme 



can be determined, where 9p is mixing angle for the states and 9yg are mixing angles for the 
decay constants /1 and /§. The theoretical and phenomenological considerations performed 
in Refs. 



16 



19 



291 ] also favor the quark-flavor basis. Therefore, we will adopt the quark-flavor 



mixing scheme to do our calculation through out the paper. 

The paper is organized as follows. In Sec. II, we outline our techniques for determining 
the 77 — 77' mixing angle, where expressions of the pseudoscalar transition form factors beyond 
the leading Fock State are provided. In Sec. Ill, we present the numerical results for the 777 
and 7/7 transition form factors, and then derive a possible range for the mixing angle by 
comparing with the present experimental data on the form factors F V1 (Q 2 ) and F^^Q 2 ). 
Some discussions of the uncertainty sources for determination are provided in Sec. IV. The 
final section is reserved for summary. 
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II. 7] - rf MIXING ANGLE AND EXPRESSIONS OF THE PSEUDOSCALAR 
TRANSITION FORM FACTORS BEYOND THE LEADING FOCK STATE 



A. Definition in the quark-flavor basis 

In the quark-flavor basis, the two orthogonal basis states are assumed to have the following 
parton composition in a Fock state description : 

, . _ \uu + dd) 

\Vq) = k + • • • > h) = *„. \ss) + ■ ■ ■ (1) 

where *f> Vi (i = q, s) denote the LC wavefunctions of the corresponding parton states, and 
the dots stand for higher Fock states. The physical meson states are related to the basis ([T]) 
by an orthogonal transformation 

/ \v ) \ f\Vq)\ /cos0 — sin0\ 

= W) • m = \ • I J • (2) 

\\T)') J \\Vs) / \sm0 COS0/ 

where <ft is the mixing angle. Under such scheme, the decay constants in the quark-flavor 



basis simply follow the pattern of state mixing due to the OZI-rule [18|, i.e. 

/ 



U(<P)dmg[f q ,f s ], (3) 



fq £s 
J rj J r) 

\ frf frf 

where the two basic decay constants f q and f s are defined as 

dx d 2 k_\ 
<Mo 1671-3 

with /i the factorization scale that is of order (9(1 GeV). 

Useful constraints to determine the 7] — rf mixing angle can be derived by considering 
the two-photon decay of rj and rf . The decay amplitudes of r\ s — > 77 and r\ q 77 have the 



/ i = 2V3/ dx ^ k J- k ± ) (4) 

J Jk*<ul 167T 3 " lV ' ±J K ' 



similar Lorentz structure as that of 7r° — > 77 30j 1 : 



^P-7l(fcl)72(te) = ^^f^lKffe)^^, (5) 

where the fine-structure constant a = 1/137, cp = (c S} c q ) = C\/2/3,5/3] for the states 

P = (rj s ,rj q ), and fp is the corresponding decay constant. Then the decay widths for 

r) — > 77 and rf —>■ 77 can be written as: 

a 2 Ml ( c q cos c s sin0\ 2 a 2 Ml, /c„sin0 c s cos0\ 2 . . 

F — = ( ~J q — md T ^-r. = W ( — + — • (6> 



1 It is noted that the higher helicity states do not have contribution to the decay amplitude. 



4 



where the two photon decay widths of rj and 77' and their masses can be found in PDG 31 ] 



and 



77^77 



Tjj'— »77 



0.46 ± 0.04 KeV, M v = 547.30 ± 0.12 MeV, 
= 4.37 ± 0.25 KeV, M v = 957.78 ± 0.14 MeV. 



From Eq.(jBj), we obtain the correlation between f q /f s and 

r_rv F„ 

U 



fs 



8vr 3 / 2 



c s a 
8^ 



\ M3 v 



COS + 



r^'— »77 



sm 1 



r^'^77 



COS< 



sm . 



(7) 



(8) 



It shows that if knowing the range of the mixing angle 0, the ranges of the decay constants 
f q and f s can be determined accordingly; and vice versa. 



B. A brief review of irj Transition form factor 



In order to calculate the 777 and 7/7 transition form factors, we first give a brief review 
of the 7T7 transition form factor F nl (Q 2 ). A comprehensive analysis of -F 7r7 (Q 2 ) has been 
given in Ref . [27J , in which the transverse-momentum dependence in both the hard scattering 
amplitude and the LC wavefunction and the contributions beyond the leading Fock state 
have been taken into consideration. Especially, a phenomenological expression to estimate 
the contributions beyond the leading Fock state has been constructed, which is based on the 
form factor's asymptotic behaviors at Q 2 — > and Q 2 — » 00. 

As has been pointed out in Ref.[23j that the transverse- momentum dependence in both 
hard-scattering amplitude and the meson wavefunction should be kept to give a consistent 
analysis of the form factor. The revised LC harmonic oscillator model as suggested in 
Ref. [321] was employed for the LC wavefunction x, k 1 ), which is constructed based on 
the Brodsky-Huang-Lepage (BHL) prescription 26j. More explicitly, the LC wavefunction 
of 7T° = -^\uu — dd) can be written as 



exp 



kj_ + m 2 



X K (m q , x,\t±) 



l(3 n 2 x(l-x) 



(9) 



with the normalization constant A n , the harmonic scale (3 n and the light quark mass m q 
to be determined. Since the contribution from the higher helicity states (Ai + A2 = ±1) 



has little contribution in comparison to the usual helicity state (Ai + A2 = 0), so we only 
write the explicit term for the usual helicity state. The spin-space wavefunction x i x i k±) 



for the usual helicity state of pion can be written as [32|, \ K ( m <?> x -> — m ql 'y m q + ^±5 
where k± = |kj_|. Furthermore, one can derive a relation between m q and j3 n by adopting 
the constraints from tc° —>■ \w and 7r° — > 77 [^tJ 

6.00^ - 1.12 + 1.3l) + 5.47 x 10 1 ) . (10) 

Jit \Pn / \ Pit J 



There are two basic type of contributions to F nl (Q 2 ) [261. l27| . i.e. 



F^m = F%\Q 2 )+Fi^m. ai) 

i^/(Q 2 ) involves the direct annihilation of (gg)-pair into two photons, which is the leading 
Fock-state contribution that dominates the large Q 2 contribution. F^ V \Q 2 ) involves the 
case of one photon coupling 'inside' the LC wavefunction of ir meson , i.e. strong interactions 
occur between the photon interactions that is related to the higher Fock states' contributions. 

By keeping the transverse-momentum dependence in both the hard scattering amplitude 
and the LC wavefunction, the valence quark state transition form factor F^\Q 2 ) can be 
written as 

Fjy\Q 2 ) = 2V3e n J o \dx] J ^^ w (x,k ± )T H (x,x' (12) 

where [dx] = dxdx'5(l—x—x f ), e n = (e 2 — e 2 ,) and the hard-scattering amplitude T H (x, x', k±) 
takes the form 

rp 1 / 1 \ q± • (z'qj. + k ± ) , 

T H {x,x,k ± ) = —3— — +{x^x). 

qi(x'q ± +k ± ) 2 

F^(Q 2 ) can be further simplified as the model wavefunction depends on kj_ through k\ 
only, i.e. ^(aCjkJ = %(x t k%_), 

V3e n r 1 dx r* 2 Q 2 
4ir 2 Jo xQ 2 Jo 

As for the second type of contribution F^ V \Q 2 ), it is difficult to be calculated in any 
Q 2 region due to its non-perturbative nature. One can construct a phenomenological model 
for Fj^ v ^(Q 2 ) based on the asymptotic behavior at Q 2 — ► and Q 2 — > 00. As suggested in 



f%\Q 2 ) = *ErL£sL M*, kl)dkl (is) 



Ref. 



271 ] , we assume it takes the following form: 



F^ V \Q 2 ) (1 + Q2/ft . 2)2 , (14) 



where k and a are two parameters that are determined by the asymptotic behaviors at 
Q 2 -> 0, i.e. 

' -■ ' " 2 " -, (15) 



a = -F^iQ )|q2_ and k 



2a 



\ -2_ 

\ dQ 2 



d rp(NV) 



7T7 



(Q 2 



where the first derivative of F^ V \Q 2 ) over Q 2 takes the form 



F(^)'(Q 2 )| Q2 _ 



"8^ 



9 /-i r 2 Q 2 /^(x,fc 2 



<9Q 2 Jo Jo 



x 2 Q 2 



From the phenomenological formula ([HP , it is easy to find that F^^Q 2 ) will be suppressed 
by 1/Q 2 to F^\Q 2 ) in the limit Q 2 -> oo. 



C. ?77 and 7/7 Transition form factors 

As for the LC wavefunctions for the pseudoscalars r\ q = ^\uu + dd) and rj s = \ss), they 
can be modeled as 13211 



.4; 



cxp 



X K (m i ,x,'k ± ) 



(16) 



where i = q, s respectively. They also depends on kj_ through kj_ only, i.e. ^ m {x, kj_) 
^(x, /c 2 ). Substituting them into the normalization (j3J), we obtain 



1 AimiPi^x{l - x) 

4 V / 2~7r 3 / 2 



Erf 



\ 8/3 2 x(l - x) 



-Erf 



\ 8ffx(l-x) 



dx 



fi 



2^3' 



(17) 



where /i stands for the factorization scale, and following the discussion in Ref . [27J , we take 
its value to be /j — 2 GeV. Under such choice, one may safely set /x — > 00 to simplify 
the computation, e.g. Erf 



,px(i-x) Imo-*oo "~ > 1) due to the fact that the contribution 
from higher |kj_| region to the wavefunction normalization drops down exponentially for the 
above model wavefunctions. 

Under the quark-flavor mixing scheme, the 777 and 7/7 transition form factors take the 
following forms: 



F m (Q 2 ) = F %7 (Q 2 )cos0-F r , s7 (Q 2 )sin, 



(18) 



and 



F V7 (Q 2 ) = F %7 (Q 2 )sin0 + F^ 7 (Q 2 )cos< 



(19) 
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where F Vql (Q 2 ) and F Vs ^(Q 2 ) stand for the r] q j and ?? s 7 form factors respectively. Similar to 
the pion-photon transition form factor, the pseudoscalar form factors Fp 7 (Q 2 ) (P = i] q and 
r] s ) can also be divided into the following two parts, 



F p ,{Q 2 )=F p V 1 \q 2 ) + F p ^>{Q 2 ). 



(NV), 



(20) 



The leading Fock-state contribution F P V J(Q 2 ) can be simplified as Eq. (fT3|) . and we only need 
to replace and ty w there to the present case of ep and \l/p, where ep = (e 2 + e 2 ,, \f2e 2 ) 
for P = (r] q , t] s ) respectively. And similar to Eq. ({T4")) . we assume the following form for the 
power suppressed non-leading Fock-state contribution F p N ) V \Q 2 ): 

a 



p ~> W) (i + g 2 /« 2 ) 2 ' 

where k and a are two parameters that are determined by: 

1 



(21) 



a 



-F Pj (Q 



^\q2^ and K 



2a 



\ dQ 2 



d W (NV) 



P 7 (Q 2 



(22) 



where the first derivative of Fp^ V \Q 2 ) over Q 2 takes the form 



? (nvy 
p 7 



(Q 2 



e P 



8tt 2 



dQ 2 Jo Jo 



5 "M, 



x 2 Q 2 



Q 2 ^0 



Naively, under strict SU(3)p symmetry, one has F Vql (Q 2 ) = F Vsl (Q 2 ), which leads to 



F V7 (Q 2 ) oc cos(0 + 45°) and F V ^(Q 2 ) oc sin(0 + 45°) 



(23) 



Therefore, if <ft varies within the region of [0°,45°] that is most probably the case, F V1 (Q 2 ) 
will decrease with the increment of 0, while F J? / 7 (Q 2 ) will increase with the increment of (ft. 
In the next section, we will show that under the case of the broken SU(3) p symmetry, such 
a fact is still exist. And then a possible range of <fi can be obtained by comparing with the 
experimental data on the transition form factors F V1 {Q 2 ) and F^^Q 2 ). 



III. NUMERICAL ANALYSIS 



With the help of the constraints from two photon decay amplitudes of ?? and??' (e^g. 
Eqs.( |71l8l) ) and the experimental data on the 777 and 7/7 transition form factors {4, 3, 6, 7|, 
one can obtain a reasonable region for <fr. There are several parameters in the wavefunction 




Q 2 [GeV 2 



FIG. 1: Color on-line: Pole-mass fit of Q 2 F V1 {Q 2 ) (Left) and Q 2 F r) ' 1 (Q 2 ) (Right) from the exper- 
imental data . The solid line stands for the average pole-mass fit with A v = 771 MeV or 
A„i = 850 MeV and the shaded band shows the experimental uncertainty. 



to be determined. As for the constitute quark masses, we take the conventional values: 
™>u, d = 300 MeV and m s = 450 MeV. By studying the pion-photon transition form factor, 
one may observe that the best fit of the experimental data is derived under the case of 



m u d ~ 300 MeV 



271 ]. For the transverse parameters j3 n , (3 q and f3 s they are proportional 



to the inverse of the charged radius of the corresponding valence quark states 33J]. The 



difference between them are less than ~ 10% as shown in Ref. 34|. For simplicity, we 
assume (3 q = (3 S = /3 n throughout this work. Uncertainties from the different choices of (3 q , 
Ps, m», d and m s will be discussed in Sec. IV. 

As for the experimental results of F^Q 2 ) and F^JQ 2 ), we take the pole-mass parameter 



fit formula that is adopted in those experiments 

1 



i, a 



QAnT [(rj/rf) -> 77] 



QVA 



2 • 

p 



As for the values of A^ and A^ 



A,, = 774 ± 11 ± 16 ± 22 MeV and A v = 859 ± 9 ± 18 ± 20 MeV 



for CLEO collaboration (4]; 



A,, = 0.70 ± 0.08 GeV and A,,, = 0.85 ± 0.07 GeV 



for TPC/Two-Gamma collaboration 



(24) 



(25) 



(26) 



Krf = 900 ± 46 ± 22 MeV 



(27) 



> 

CD 
O 

o 



o 




Q 2 [GeV 2 



FIG. 2: Color on-line: Q 2 F rjl {Q 2 ) and Q 2 F r] i 1 (Q 2 ) with BHL-like wavefunction. The dash-dot line, 
the solid line and the dashed line are for = 37.0°, cf> = 38.0° and eft = 39.0° respectively. It shows 
that Q 2 F V ry(Q 2 ) decreases and Q 2 F v i 1 (Q 2 ) increases with the increment of cj>. The shaded band is 
the region allowed by the experiments % 



mcrea, 1 

m. 



for L3 collaboration 



and 



A,, = 0.84 ± 0.06 GeV and A v = 0.79 ± 0.04 GeV 



(28) 



for CELLO collaboration jlj]. Averaging the above experimental values, we obtain the center 
value for A^ and A v >, i.e. A,, = 771 MeV and A v > = 850 MeV. We draw the pole-mass fit of 
the form factors Q 2 F m (Q 2 ) and Q 2 F rj i 1 {Q 2 ) in Fig. ([I]), where the shaded band is derived by 
adopting the pole-mass fit formula fl24l) and by varying A,, and A v i within the widest possible 
range allowed by the above experimental results 2 . The shaded band (region) for Q 2 F rjl {Q 2 ) 
and Q 2 F V ^(Q 2 ) can be regarded as constraints to determine the rj/r)'— wavefunctions, i.e. 
the values of the parameters in the wavefunctions and also the mixing angle <fr should make 
Q 2 F m (Q 2 ) and Q 2 F V ^(Q 2 ) within the region of the shaded bands as shown in Fig. (CD). 
From Eq. fll7p . we obtain 



A q ~ 2.77 x 10 2 f q and A s ~ 2.85 x 10 2 / s 



and then with the help of Eqs. ( 1711811181)191) . it can be found that only the mixing angle <ft is 
undetermined. As shown in Fig.(j2J), one may observe that the value of Q 2 F m (Q 2 ) decreases 



2 There we do not take the weighted average of these experiments and treat them on equal footing, as these 
experiments are concentrate on different energy regions and only few data are available. 
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with the increment of 0, while the value of Q 2 F v r 1 (Q 2 ) increases with the increment of 4>, 
and such a fact can be used to determine a possible range for <fi by comparing with the 
experimental data. In fact, it can be found that the lower limit of <fi is determined by 
Q 2 F rjl (Q 2 ) and the upper limit of <fi is determined by Q 2 F V ' 1 (Q 2 ), i.e. 

(p ^ 38.0° ± 1.0°. (29) 

Furthermore, we obtain j- = 1.07 ± 0.01 and j- = 1.24 ± 0.10. By using the correlation 



18 



281 ]. we obtain 



between the quark-flavor mixing scheme and the octet-singlet scheme 

Q p = 0-arctan v / 2 = -17.0° ± 1.0° , (30) 
e 8 = 6- arctan 4^ = -20.7° ± 1.0° 1 B 1 = 6- arctan 4^ = -12.6° ± 3.3° (31) 

Jq Is 

and 



A = yJIZ+lIL/^ = L19 ± o 07 ) A = ^ 2 J1±H/U = 1.13 ± 0.03 , (32) 



IV. DISCUSSION ON THE UNCERTAINTIES OF DETERMINING <p 

In this section, we discuss the uncertainties for determining the mixing angle from the 
above approach. First, we compare the differences caused by the different model wavefunc- 
tions, e.g. the BHL-like one and the CZ (Chernyak-Zhitnitsky) -like one 3s| which is much 
broad. And then, we restrict ourself to use the BHL-like model wavefunction for a detail 
analysis on the effects to <ft determination caused by each uncertainty sources separately, 
where the uncertainty sources mainly include the value of f3 q and f3 s , the masses of the con- 
stitute quarks u/d and s, and the possible intrinsic charm components in rj and rj'. Some 
other even smaller uncertainty sources for the electro- magnetic transition form factors such 
as the gluon component in rj/rj' will not be discussed 3 . 



A. Model dependence 

One typical broad wavefunction is described by the CZ-like wavefunction. For convenience 
and simplicity, we take m q = 0.30 GeV and m s = 0.45 GeV, and (3 q = (3 S = j3 n . It can be 



3 The gluon contributions might be important to some other exclusive processes like B meson two-body 
non-leptonic exclusive decays. And a discussion on the two gluon components in the form factors can be 
found in Ref . [361] . 
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found that f3 n = 0.70 GeV for the case of CZ-like wavefunction 27]. 
The CZ-like wavefunctions for r) q and r) s take the form 



^f(x,k ± ) = A? z (l-2x) 
where i = q, s. Following a similar procedure, it can be found that 



/ k 2 ± +m 2 \ K 
exp k— v mj, x, k i 



(33) 



~ 8.86 x 10 2 / 9 and A c s z ~ 8.95 x 10 2 / s • 

Through numerical calculation, it can be found that the 777 and 7/7 transition form factors 
from the CZ-like wavefunction raise faster than the case of the BHL-like wavefunction. And 
the shapes of the 777 and 7/7 form factors from the BHL-like wavefunction are more close to 
the pole-mass parameter fits of the experimental data. Although the model is different, one 
may find that the value of Q 2 F m (Q 2 ) decreases with the increment of 0, while the value of 
Q 2 F r/ i 1 (Q 2 ) increases with the increment of 0. So the range for can also be estimated, i.e. 
= 38.0° ± 1.0°, which is close to the case of BHL-like wavefunction. This shows that the 
77 — 77' mixing angle is almost model-independent in our approach. 

Recently, BaBar collaboration has measured the value of Q 2 F m (Q 2 ) and Q 2 F ri i 1 (Q 2 ) at 



Q 2 = 112 GeV 2 [lOfl: Q 2 F rn (Q 2 ) = 0.229 ± 0.030 ± 0.008 GeV and Q 2 F V7 (Q 2 ) = 0.251 ± 



O 2 F 1 (O 2 ) , 

0.019 ±0.008 GeV, and the ratio of the form factors k = q2^ ( q 2) \q2 =112 GeV 2 = 1-10 ±0.17. 
Under the case of G [37°, 39°], for the BHL-like wavefunction, we have Q 2 F m (Q 2 ) = 
[0.176, 0.190] GeV and Q 2 F V ^(Q 2 ) = [0.228, 0.277] GeV at Q 2 = 112 GeV 2 , which is close to 
the experimental values 4 and leads to k — 1.44 ± 0.06. While for the CZ-like wavefunction, 
we have Q 2 F m (Q 2 ) = [0.267,0.297] GeV and Q 2 F^ 7 (Q 2 ) = [0.388,0.411] GeV at Q 2 = 
112 GeV 2 , which is somewhat bigger than the experimental values and leads to k — 1.41 ± 
0.06. As a comparison, one may conclude that the asymptotic behavior of the form factors 
Q 2 F rjl (Q 2 ) and Q 2 F r) > 1 (Q 2 ) disfavor the CZ-like wavefunction but favor the asymptotic like 
wavefunction. 

It is due to these differences that the form factors of CZ-like wavefunction and the BHL- 
like one are affected differently by the following considered uncertainty sources. For example, 



as shown in Ref. 27J], the best fit of 717 form factor to the experimental data is obtained with 



4 



It is also close to the theoretical predictions based on the asymptotic wavefunction that has been clearly 
shown in Fig. 13 of Ref.0. 
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m q ~ 300 MeV in the case of BHL-like wavefunction, which is shifted to m q ~ 400 MeV for 
the case of CZ-like wavefunction. 

In the following, we will only take the BHL-like wavefunction as an explicit example to 
show the uncertainties and the case of CZ-like one can be done in the similar way. For clarity, 
in studying of the uncertainty caused by a certain source, the other uncertainty sources are 
taken to be their center values as adopted above. 



B. Uncertainty A(j) m from m q and m s 

We take a wider range for m q and m s to study their effects to the mixing angle 0, 
e.g. m q = 0.30 ± 0.10 GeV and m s = 0.45 ± 0.10 GeV. Under the present case we adopt 
j3 p = (3 S = j3 n , where the value of @ n varies within the region of [0.48, 0.70] GeV according to 



the value of m q 



271 ] . From Eq. (jl7p . we obtain the uncertainty from the constituent quark 



masses 

7+1.00 v> -i r\2 £ _J A ^ o OK+1-99 w 1fl2. 



A p ~ 2.77™ x 10V P and A s ~ 2.85™ x 107.. 



S ) 



where both A p and A s increase with the increment of m q and m s respectively. And, it can 
be found numerically that 

A0 m < ±0.5°. (34) 



C. Uncertainty Ac/r from (3 q and f3 s 

Due to SU(3)-symmetry breaking, there are differences among (3 q and (3 S , which is 
smaller than 10% through a light-cone quark model analysis 34|. For clarity, we choose 
broader ranges (3 q = 0.55 ±0.10 GeV and (3 S = 0.55 ±0.10 GeV to make a discussion on how 
these transverse size parameters affect the mixing angle. Other wavefunction parameters 
are fixed by setting m q = 0.30 GeV (or equivalently @ n = 0.55 GeV) and m s = 0.45 GeV. 
Under such condition, we have the uncertainty from the transverse parameters (3 q and /3 S , 

A p ~ 2.77^;^ x 10 2 / p and A s ~ 2.85^;^ x 10 2 / s , 

where both A p and A s increase with the decrement of (3 p and /3 a , respectively. And, it can 
be found numerically that 

Acp 13 < ±2.0°. (35) 
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D. Uncertainty A0 C from the intrinsic charm component 



It has been suggested that a larger intrinsic charm component might be possible to explain 
the abnormally large production of rf in the standard model 37] . However, some studies in 



the literature disfavor such a large portion of intrinsic charm component, e.g. [16l.l21l.l38l.l39] 
and references therein. It is has been found [3] that the mixing between the cc state with 
qq-ss basis is quite small, e.g. less than 2%. So for simplicity, we do not consider the mixing 
between cc and qq-ss basis. And then, we have 

F m (Q 2 ) = F Vql (Q 2 ) cos0 - F„ s7 (Q 2 ) sin0 + F^Q 2 ) (36) 
F^Q 2 ) = F^Q 2 ) sin0 + F VsJ (Q 2 ) cos0 + F^'jQ 2 ), (37) 

where F^(Q 2 ) and F^(Q 2 ) corresponds to the contributions from the intrinsic charm 
component in r\ and rj' respectively, which will be calculated in the following. 

The wavefunction of the "intrinsic" charm component rj c = |cc) can be modeled as 

k 2 ; + m 2 



A 



v/v' 



exp 



l(3 2 x(l -x)_ 



X (m C} x,k ± ) 



(38) 



where we adopt (3 C = {3 n = 0.55 GeV and m c = 1.5 GeV 5 . The overall factor A° is 
determined by the wavefunction normalization similar to Eq.t fiTl) . which shows 



A c v = 9.45 x 10 3 ^ and Afy = 9.45 x 10% c , , 



3 t c 



where /f and /£> are related through [16] 



f c 



arctan 



V2f s 



(39) 



and if taking <fi = 38.0° ± 1.0°, we have 



ft 



[0.36,0.40]. 



It is noted that different from the 7^7 and 7/ s 7 transition form factors, the helicity-flip 
amplitude that is proportional to the current quark mass cannot be ignored for the present 
case. For the rjcj transition form factor, a direct calculation shows 6 



2V6e 2 







[dx] 



d 2 k, 



16tt 3 



(40) 



5 By varying /3 C and m c within their possible regions, the following results will be slightly changed and the 
present conceptional results are the same. 

6 There we will not consider the non-valence charm quark states' contribution since it is quite small due to 
the large charm mass effect. 
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FIG. 3: Q 2 F m (Q 2 ) and Q 2 F V ^{Q 2 ) 
shows the experimental uncertainty 4 



p = 38.0° with BHL-like wavefunction. The shaded band 
I, 0, Q]. The dash-dot line, the solid line, the dashed line 



and the triangle line are for fZ = 0, —5 MeV, —15 MeV and —50 MeV respectively. 



where the hard-scattering amplitude Tw(x, x',kj_) that includes all the helicity states (Ai + 
A 2 = , ±1) of rj c takes the form |2ll, kdj. 

q_L ■ (x'q_L + kj.) 



nA X : X J k±) — r . , , 

q]_[(x'q_L + k_l)^ + m z c \ 



+ (X <-»• X ). 



(41) 



The above formula can be further simplified by doing the integration over the azimuth angle: 



1 + 



1- z-y 2 



(z + (1 - y) 2 )(z + (1 + y) 2 ) 



h i dk i 



(42) 



where z 



x 2 Q 



% and y = ^. 



Taking = 38.0°, we draw in Fig. ([3]) how the value of f5 affects the form factors 
Q 2 F V1 (Q 2 ) and Q 2 F V >^(Q 2 ). One may observe that the experimental data disfavors a larger 
portion of charm component as \f^,\ > 50 MeV. Or inversely, it can be founded that under the 
condition of \f^,\ < 50 MeV, the uncertainty from the possible intrinsic charm components 
is given by 

A0 C < ±1°. (43) 



E. A summary remark on (ft 



Under the quark-flavor mixing scheme, and by carefully dealing with the higher Fock 
states' contributions, a possible range for can be derived by comparing with the exper- 
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imental data on the form factors F m (Q 2 ) and F r) i 1 {Q 2 ) due to the fact that the value of 
Q 2 F V1 {Q 2 ) decreases, while the value of Q 2 F v t 1 (Q 2 ) increases, with the increment of in 
the whole Q 2 region. It has found that the allowable range for the mixing angle is 

0^38.0° ±1.0° ±2.0°, (44) 

where the first error is from experimental uncertainty and the second error is from 

the uncertainties of the wavefunction parameters and the possible intrinsic charm component 
in 7] and rj . It should be note that the second uncertainty is lower than the direct sum of 
the errors caused by each source separately, i.e. A0 m + A0^ + A0 C , which is due to the fact 
that these uncertainty sources are correlated to each other. 

V. SUMMARY 

In the present paper, we have performed a light-cone pQCD analysis of the 777 and 7/7 
transition form factors F m (Q 2 ) and F r/ / 7 (Q 2 ) involving the transverse momentum corrections, 
in which the rj — rj mixing effects and the contributions beyond the leading Fock state have 
been taken into consideration. For such purpose, we have adopted the quark-flavor mixing 
scheme for the 77 and rj mixing and have constructed a phenomenological expression to 
estimate the contributions beyond the leading Fock state based on its asymptotic behavior 
at Q 2 — > and Q 2 — ► 00. It has been found that the value of Q 2 F V1 (Q 2 ) decreases, while 
the value of Q 2 F r) i 1 {Q 2 ) increases, with the increment of in the whole Q 2 region, and then 
a possible range for can be determined by comparing with the experimental data, which is 
= 38.0° ± 1.0° ± 2.0° with the first error coming from experimental uncertainty |4j, |5|, |6j, |7] 
and the second error coming from the uncertainties of the wavefunction parameters and the 
possible intrinsic charm component in 7/ and rj . A more accurate weighted average of the 



values as described in 
and = 38.8° ± 1.2° 



above mentioned value together with the seven adopted experimenta 

Jlef. [1^1, and the two new experimental values, = 41.2° ±1.2° [l2 

13j . yields = 39.5° ± 0.5°. Furthermore, our results show that the 7/ — rj mixing angle 
depends on the different wavefunction models slightly. However the asymptotic behavior of 
the form factors Q 2 F m (Q 2 ) and Q 2 F V ' 1 (Q 2 ) disfavor the CZ-like wavefunction but favor the 
asymptotic like wavefunction. It has been found that the intrinsic charm component in 77 and 
rj can not be too big, e.g. \fZ\ < 50 MeV. Such a conclusion agrees with other investigations 
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l2ll . l38l . |39j. These results are helpful to understand other exclusive processes involving 
the pseudo-scales r\ and rj '. 
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